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Abstract 

The holographic model for S-wave high T^ superconductors developed by Hartnoll, Herzog 
and Horowitz is generalised to describe D-wave superconductors. The 3+1 dimensional grav- 
itational theory consists of a symmetric, traceless second-rank tensor field and a U{1) gauge 
field in the background of the AdS black hole. Below Tc the tensor field which carries the U{1) 
charge undergoes the Higgs mechanism and breaks the U{1) symmetry of the boundary theory 
spontaneously. The phase transition characterised by the D-wave condensate is second order 
with the mean field critical exponent /3 = 1/2. As expected, the AC conductivity is isotropic 
below Tc and the system becomes superconducting in the DC limit but has no hard gap. 



I. INTRODUCTION 



One of the unsolved mysteries in modern condensed matter physics is the mechanism 
of the high temperature superconducting (HTSC) cuprates These materials are lay- 
ered compounds with copper-oxygen planes and are doped Mott insulators with strong 
electronic correlations. The pairing symmetry is unconventional and there is a strong 
experimental evidence showing that it is D- wave [21. It is speculated that the pairing be- 
tween electrons is mediated via strong anti-ferromagnetic spin fluctuations in the system. 
A prominent strong coupling theory is proposed by Anderson, called the resonant valence 
bond (RVB) theory, which describes liquid state with spin- singlets. Upon hole doping, the 
Neel order is destroyed and give rise to sup er conduct ivityjsl. Several gauge theories have 
been proposed to formulate the RVB physics, by enforcing the double occupation con- 
straint in the strong coupling limit jj]. The problem is difficult due to the strong-coupling 
nature of the theory. Although significant progress has been made in the past few years, 
alternative approaches may be valuable to tackle the problem. 

One alternative approach is the holographic correspondence between a gravitational 
theory and a quantum field theory, which first emerged under AdS/CFT correspondence 
SKZl . This method has provided a useful and complimentary framework to describe strong 
interaction systems without a sign problem (see e.g. js iSl). In the original top-down ap- 
proach, both the gravity side and the field theory side of the theories are precisely known. 
Later applications assume that the correspondence exists among different pair of theories 
and try to make predictions from one side of the correspondence. More specifically, in this 
bottom-up approach, usually the gravity side of the theory is explicitly constructed with 
the desired symmetries, then physical observables (matrix elements) of the field theory 
side are predicted through the above mentioned correspondence. 

Recently, a grav itational model of hairy black holes jiol, have been used to model 



S-wave HTSC 



18 



21| . In those class of models the Abelian symmetry of a complex 



scalar field is spontaneously broken (i.e. the Higgs mechanism) below some critical tem- 



perature. 



he Meissner effect was soon observed by including magnetic field in the 



background 22, l23[ . The effect of superconducting condensate on the holographic fermi 



surface has been studied by calculating fermionic spectral function 24l-l26| . Interestingly 



the properties of spectral function appeared to have similar behaviour to that found in the 
angle resolved photo-emission experiment. Motivated by all of these s-wave studies, holo- 
graphic dual to P-wave superconductor has been proposed by coupling a SU (2) Yang- Mil 



field to the black hole, where a vector hair develops in the superconducting phase 27H30|. 



Behaviour of fermionic spectral function has also been studied in those p-wave supercon- 
ducting background 3l| . So far, the bottom- up construction of holographic superconduc- 
tor assuming the existence of gauge/gravity duality has been discussed. However, in the 
string theory framework people also have studied top-down _approach considering various 
D-brane configurations in the AdS black hole background 32 1. 

In this work, we try to construct a minimal grav itational model that models D-wave 
HTSC We replace the complex scalar field in 18| by a tensor field whose condensate 
breaks the symmetry spontaneously below Tc and the condensate becomes zero and the 
symmetry is restored above T^. The critical exponent /3 gives the mean field value 1/2. 
The real part of the conductivity computed from linear response has a delta function 
at zero frequency which corresponds to static superconductivity below Tc. Above Tc, 
the delta function disappears as expected and the conductivity becomes constant in fre- 
quency. It is expected that there is no "hard gap" in the real part of the conductivity 
and the conductivity should be isotropic even though the condensate is not (for a model 



calculation, see 



33|). Both features are seen in our results. 



II. A HOLOGRAPHIC MODEL FOR D-WAVE HTSC 

Our goal is to consider a minimal (3+1 dimensional) holographic model that gives 
rise to (2+1 dimensional) D-wave superconductivity. The construction will be similar 
to that of the S-wave case 13] with a spontaneous local f/(l) symmetry breaking in 
the bulk leading to a spontaneous breaking of global U{1) symmetry at the boundary. 
Thus, strictly speaking, the boundary theory is a super-fluid. One can still study the 
current-current correlator which could be interpreted as the conductivity. 

To have a D-wave condensate at the boundary, we introduce a charged tensor field 
in our dual gravity theory. Assuming the D-wave condensate originating from electron- 
electron pairing, the D-wave nature gives a symmetric wave function for the pair. So, wave 
function for this electron pair has to be a spin singlet such that its total wave function 
is anti-symmetric. A 3 x 3 symmetric traceless tensor has 5 components which can be 



used to describe a D-wave state. We will promote this symmetric traceless tensor field 



to include time components and denote the field as B^,y {fi, u = 0,1, 2, 3), i.e. B^^ = B, 



UfJ, 



and B^ = 0. However, it is important to note that the interacting higher spin fields, in 
general, require to satisfy additional constraints in addition to the equations of motion 
to remove the unphysical degrees of freedom, see the discussion in 3J-|36| for example. 
Observing that there is no available consistent model for our purpose, we would like to 
propose a truncated model which has sufficient ingredients to catch some features of D- 
wave superconductor. It would be an important but difficult task to construct a complete 
theory which we would like to postpone for the future in order to attain our simple goal 
through the present exercise. 

The desired complete action including gravity, U{1) gauge field, tensor field and other 
auxiliary fields, may take the following formjs^]: 

6 



S 



r = 



1 

2 



{D^B^^yO^B 



fiTJU-y , _2 



4 1"^-^ 



where R is the Ricci scalar, the term gives a negative cosmological constant and L 
is the AdS radius which will be set to unity in the units that we use. k,^ = SttGn is the 
gravitational coupling. is the covariant derivative in the black hole background {D^ = 
dn + iAfj^ in fiat space), and q and are the charge and mass squared of i?^^, respectively. 
The terms associated with auxiliary fields are included in Ca- This Lagrangian is the same 
as that appears in 34] in the fiat space limit. However in the subsequent studies after 34] 
it appeared that construction of higher spin field Lagrangian coupled to gravity or f/(l) 



gauge in a gauge invariant way is non-trivial even in Minkowski space [38[. But later on 
in the context of gravitational interaction 39| it has been shown that this task may have 
a solution in AdS space. Also there exists a recent attempt to construct gauge invariant 
f/(l) charged massive spin 2 particles in AdS space at linear approximation 40|]. In spite 
of all these studies fully consistent formulation of interacting higher spin gauge fields is 
still lacking. Keeping this in mind that the action in ([1]) without Ca contains spurious 
degrees of freedom that may cause the instability. We will proceed by assuming that the 
constrains can eliminate the instabilities but still allow the D-wave condensation. 

Cm might look more familiar with the rescaling B^^ — > g-B^j^ and — )■ qA^. Here we 



also concentrate on the "probe limit" 18|] where the back-reaction to the background can 



be ignored. This limit is exact when g — )■ oo. In the probe hmit, can be treated as a 
perturbation on top of the 3+1 dimensional AdS black hole background: 

ds'^ = -g(r)dt^ + —- +r^(dx^ + dy^), (2) 
9[r) 

3 

where g{r) = r'^ — ^ and tq is the horizon size. The Hawking temperature for this black 
hole T = ^. 



47r 

As in S-wave case [18|], electric field can exist in the bulk by the appropriate choice 
of boundary conditions. The charged tensor field, which can be considered as charged 
particles, experiences a force under the electric field, with positive (negative) charges re- 
pelled(attracted) away from(toward) the black hole. One the other hand, the black hole 
tries to pull all the charged particles in it. At lower T, the black hole is smaller and the 
gravitational pull is weaker. Thus, the positively charged particles have a bigger chance 
to stay outside the horizon and form the condensate. At very large T, the gravitational 
force from the large black hole is strong enough to pull all the charged particles into the 
horizon such that there is no condensate. Thus, we have a phase transition. 

We are interested in describing the D-wave SC in the continuum such that there is a 
condensate on the x-y plane on the boundary with translational invariance. Rotational 
symmetry is broken down to Z{2) with the condensate changing its sign under a 7r/2 
rotation on the x-y plane. To incorporate these features, we use an ansatz for the Bfj,^, 
and the gauge field A^, i.e., 

= diagonal (0,0, /(r),-/(r)), A = (j){r)dt. (3) 

After plugging in this ansatz, we have the equation of motion for B 

L 9{r) J I 9{r)^ g{r) } 

and the corresponding Maxwell's equation is 

2 ,/// N , o ^ 4/^(r)0(r) 
r (r) + 2r0 (r) 5—- — = 0, (5) 

where the ' is the derivative with respect to r. 

We would like to choose the solution such that 0(r) has the asymptotic form 

0(r)^/i+- (6) 
r 



near the boundary (r — oo), where /i is interpreted as the chemical potential and p as 
the charge density in the boundary theory. Here, we will first assume this and then show 
that indeed this can be satisfied later. Now, the asymptotic of EqJl]) has the asymptotic 
form 



rV"(r) - (m2 + 4) /(r) ^ (7) 
near the boundary, which yields 

A, = l±^E±^. (8) 

If we interpret /o as the source and /i as the vacuum expectation value (VEV) of the 
operator that couples to B at boundary theory, we need m? > —4 (and A_ < 0) such that 
the /i term is constant or vanishing at the boundary. After setting the source /o = and 
using A_ < 0, Eq.Q indeed gives the asymptotic solution of Eq.([6]). Note that the /o?"^+ 
term does not impose a constraint on by requiring that the third term on the LHS of 
Eq.Q to be smaller than the other two terms since we have imposed fo = 0. One way 
to see this is to do the integration of the differential equations from the boundary, then 
/(r) — > /ir^", 0(r) fi + ^ satisfy the asymptotic behaviors of Eqs.(jl]) and 1^. The 
order parameter of the boundary theory can be read off from the asymptotic behavior of 

m = ( I (9) 

where are the indexes in the boundary coordinates {x,y). 

It is useful to note that the action and the equations of motion are invariant under the 
scaling 

(t,r,x,y) {t/c,cr,x/c,y/c) , 
(ro,T,5((r)) (cro, cT, c^5((r)) , 

(/(r),0(r)) -> {c'fir),c<P{r)). (10) 

Thus, we can always scale /i — 1. This also helps to keep track of the scaling dimension for 
observables, e.g., the scaling dimensions for /i, p, and fi are 1, 2, and 2 — A_, respectively. 
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FIG. 1: (color online) The dimensionless D-wave condensate f\l shown as a function of 
T jTc for 771^ = —1/4. The condensate goes to zero at T = Tc oc /x. The critical exponent 
/i — )• c [Tc — T)^ for Tc — T — )• 0^ is of the mean field value /3 = 1/2. 

In Fig. 1, we show the numerical result between the dimensionless quantities /i/ jJ^"^^ 
and T/Tc for m? = —1/4. It is a second order phase transition. Numerically, the critical 
exponent (3 defined as /i — )■ c (T^ — T)'^ for Tc — T — 0^ is very close to the mean 
field value /3 = 1/2. Below we show that only the values /3 = 1/2,3/2,5/2,... satisfy 
the equations of motion. Thus, without fine tuning, one would get the mean field value 
(3 = 1/2. 

Now we present the derivation. The metric g{r) is a smooth function of e = — T, 
while and / can be expanded as 



g{r,T) = agir) + bg{r)e + 0{e^), 
0(r,T) = e'{a4r) + b^ir)e + 0{e')), 

/(r,T) = e-{afir) + bf{r)e + 0{e')). (11) 

Since, Eq.(jlj) is a linear equation in / and the pre- factors of f"{r) and f'{r) are polynomials 
of e, the pre-factor of /(r) has to be a polynomial of e as well in order the satisfy the 
equation. This implies k is an integer. At and at the boundary, gives the value of 
chemical potential which is finite. This yields k = 0. 

Analogously, the Maxwell equation, Eq.([5]), is linear in (p. The pre-factors of (f)"{r) and 
(/)'(r) are polynomials of e and thus the pre-factor of 0(r) is required to be a polynomial 
of e. This yields 2n to be an integer. We also know that n > for a second order phase 



transition, and hence only /3 = 1/2, 3/2, . . .are allowed. 



III. CONDUCTIVITY 

In this section, we compute the conductivity of this D-wave HTSC by linear response. 
The conductivity tensor aij can be defined through the linear response relation 



(T- -E- 



(12) 



where z,j = 1,2,. J and E are the electric current and electric field, respectively. Fol- 
lowing the approach of 



we perturb the gauge field by 6A = e '^^^Ax{r)dx. To get 



a consistent set of equations, we also need to perturbe 6Brx = SBj-r = iKxi^)^ and 
xt = btx{'r)e respectively. The resulting equations of motion are 



gA'i + g'A'^ + 



4/^ 



A. = 0, 



(13) 



gbrx" + 2g'brx' 



3__ 3__^9_ 



{CO - 4>y 



2/ 



Ax 



9' 



Hx, 



gbtx" + 



g" g' g {u - 0)- 

2 r g 



btx = (w - (l))g'br 



(14) 



(15) 



In principle, we can also add other 5B^i, components in the perturbation. However, 
those components do not couple to 5A^ to the quadratic order in the action. So, if we 
set the initial condition of our system to be in the ground state before 5^4^ perturbation 
being turned on then those extra B field perturbations will not be produced. However, 
in the full stability analysis, those 5A^ independent perturbations are important. We will 
defer this stability analysis to future study. 

Eq.(|T3l) is very similar to the S-wave case and is decoupled from 5-6^^. Near the 
boundary, we have 



r^A" + 2rA'^ ~ 0, 



which yields the asymptotic form 



Ax Ax.Q + 



AxA 



(16) 



(17) 



where A^^ is the x-component gauge field at the boundary whose time derivative gives 
E^, and A^^i is the expectation value of the current operator J^. The ratio of and 
is the frequency dependent conductivity 

cr (u) = a^^ (u) = _l^pL_ (18) 
The fact that Eq. (fT3l) depends only on implies 

(Ty^ (oo) = 0. (19) 

This is dictated by the refiection symmetry with respect to the y = plane. 

We are now focusing on the case > —2, where the asymptotic forms of Eqs. fll4l) 
and (fT5i) are particularly simple: 



brJ' + ArbJ - (m^ + 6) Kx ^ 0, (20) 
- (m^ + 4) ^ 0. (21) 



These two equations can be solved with 

brx hrxfiT^^ + Kx^T^' 

btx ^ btx,or'^+ + btx,ir'^- , (22) 

where A± = -3±V33+4 m^ defined in Eq.®. Here, we also identify brx,o and btx,o 

as the source terms and brx,i and btx,i are the normalizable fiuctuations. 

Near the horizon, g{r) = Sr^dr + 0{dr'^) with dr = r — Tq and (j){r) = 0{dr). The 
equations of motion become 

9dr^A'^ + 9drA'^ + -^Ax = 0, (23) 



r, 



9dr%rx" + ISdrbrx' + ^brx = -^hx " ^^^dvAx, (24) 

Tq r^dr 

9dr\x" + —btx = 9ujdrbrx- (25) 



The solutions near the horizon are 



Ax ax,idr '^-o j^- a^^^dr^'-o , 

btx 3iro(brx^dr~'^^^ -brx,2dr~'^ -brx,3dr^>^^ + KxAdr^^ , 

brx brx,idr ^'■0 + brx,2dr + brx,3dr '-"-o + brx,4dr s-o . (26) 



Refer] Re[cr],Im[cr] 

1-5 r 




FIG. 2: (color online) The real (left plot) and imaginary (right plot) of conductivity shown as a 
function for frequency oj for different temperatures. Above Tc, Re[a (w)] = 1 while Im[a (a;)] = 0. 
Below Tc, Re[(j (oj)] has a 5 (oj) delta function whose height decrease in T and vanishes at Tc- 
The right most curve has the lowest T, which implies the zero temperature gap cOg/Tc — 13. 



(The construction in 



lg | for S-wave gives the value 8 for this gap.) 



The ingoing wave boundary condition 



41 



42], which sets the wave falling into the 



horizon, demands ax,2 = brx,3 = Kxa = 0. We further set the divergent term brx,2 = to 
keep the action finite. Now, 6ra;,o{i) and btx,Q{i) in the Eq.( l22|) are linear combinations of 
0^,1 and brx,i- So we have both normalizable and non-normalizable solutions for btx and 
brx perturbations. The divergent source terms b^x.o and btx^ near the boundary can be 



cancelled by counter terms 



19| 



In Fig. 2, we plot the real and imaginary part of a (u) for different T. The behaviours 
are similar to that of S-wave HTSC Re[a (u)] has a delta function behaviour at w = 
corresponding to infinite DC conductivity when T < Tc. On the othe hand for T > Tc, 
the delta function and Jm[(j (u)] disappear and Re[a (u)] becomes u independent. There 
is no "hard gap " in our dual boundary superconducting system because Re[a {u)] does 
not vanish even for arbitrary small u. One can also read off this soft gap from the plot. 



i.e. oOg/Tc ~ 13. It is larger than the one obtained in the construction for S-wave |l SI], 
where oOg/Tc ~ 8. This may imply our D-wave pairing requires higher energy than the 
S-wave one. 

Unlike the case for S-wave superconductor, the vanishing of the gap is actually expected 



in the D-wave case. In the BCS-type theory (see, e.g. 43|), the lowest dimensional D- 

wave operator for two fermion pairing is Oij = tp'^ ( d i d j — d '^5ij/2] ip, where d « = 
— )■ ^ 

d i — d i is the relative momentum between the two fermion which is invariant under 



Galilean transformation. The leading order Lagrangian in a weakly interaction theory is 

£ = £0 - c (0{, + 4) (O,, + J,.j) + cJj, (27) 

where Cq is the free Lagrangian, c is the coupling and Jij is an external source. Under a 
Hubbard Stratanovich transformation, the Lagrangian can be rewritten as 



where B and B* are auxiliary fields. After integrating over the auxiliary fields, C is 
recovered from C. It is clear that (Oij) in C is {Bij) in The gap equation of C gives 
the dispersion relation 



^'^ = \|{^-^') +\B^Jhk/■ (29) 

The gap \Bijkikj\ oc \kl — ky\ vanishes at four nodes kl = ky. So, naturally gapless 
excitations can contribute to conductivity. This makes the conductivity for a D-wave 
superconductor gapless. In the S-wave case, however, the gap is isotropic and does not 
vanish in any direction leading to a hard gap in conductivity. 

If we change the gauge field perturbation to 6A = e"""^* {A^{r)dx + Ay{r)dy), then 
there will be response from 6Brx,SBt;j.,6Bry and 6Bty. A{r) = Ax{r)x + Ay{r)y satisfies 
the same differential equation as Eq. (fT3|) : 

9A" + g'A'+(^^-^^A = 0. (30) 

This shows that the conductivity is isotropic: 

aij{u) = (r{u)Sij. (31) 

This might seem surprising at the first sight because the condensate is not isotropic. 
However, this is a consequence of the symmetries that o".y has in the D-wave case. In 
the linear response theory, aij is a current-current correlator which can be schemati- 
cally denoted as aij ~ (fi |[Jj, Jj]| fi), where the matrix element denotes an ensemble 
average. Under a 7r/2 rotation along the z-axis [R), R~^JiR = ^ijJj-, where ejj is an 
anti-symmetric tensor, and assuming the ensemble average is governed by properties of 
the ground state which has the condensate structure of Eq.(l9]), so -R jfi) = — 1^2). Then, 
aij ~ (fi |[Ji, Jj\\VL) = {Q \R~^ [Ji, Jj] R\ Q) = {Q \ [eikJk, ejiJi]\Q). This implies a^x = (^yy 
and axy = —ayx- A similar analysis with parity operator with respect to the x-axis gives 



o'xy = <^yx = 0. Thus, we have aij oc 5ij. An explicit microscopic model calculation 
also yields an isotropic conductivity for a D-wave superconductor. 



33| 



IV. CONCLUSION 



We have constructed a minimal holographic model for high Tc D-wave superconductors. 
We follow closely the work of Hartnoll, Herzog and Horowitz on the S-wave case. The 
3+1 dimensional gravitational theory consists a symmetric, traceless second-rank tensor 
field and a U{1) gauge field in the background of the AdS black hole. Below the 
tensor field is Higgsed to break the f/(l) symmetry at the boundary theory. The phase 
transition characterised by the D-wave condensate is second order with the mean field 
critical exponent /3 = 1/2. As expected, the AC conductivity is isotropic; below Tc, the 
system becomes superconducting in the DC limit but has no hard gap. 
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